UWThPh-2001-15 
April 2001 



Isospin Violation 
and the Magnetic Moment of the Muon^ 

V. Cirigliano, G. Ecker and H. Neufeld 



Institut fur Theoretische Physik, Universitat Wien 
Boltzmanngasse 5, A-1090 Vienna, Austria 



Abstract 

We calculate the leading isospin-violating and electromagnetic corrections for the decay 
t~ — > 7t 7t~v t at low energies. The corrections are small but relevant for the inclusion 
of r decay data in the determination of hadronic vacuum polarization especially for the 
anomalous magnetic moment of the muon. We show that part of the systematic differences 
between the measured form factors in t~ — > 7i 7T^u T and e + e~ — > 7r + 7r~ is due to isospin 
violation. 
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1. Precise knowledge of hadronic vacuum polarization is essential for a reliable deter- 
mination of both the running of the QED fine structure constant and of the anomalous 
magnetic moment of the muon a M . For the latter, the low-energy structure of hadronic 
vacuum polarization is especially important. In fact, about 70 % of a™ cpo1 is due to the 
two-pion intermediate state for AM 2 < ^ < 0.8 GeV 2 (see, e.g., Ref. 0). 

A precision of 1 % has been achieved in the calculation of a™ cpo1 by including [0] the 
more accurate r decay data [[J in addition to experimental results for <r(e + e~ — > hadrons). 
This is possible because of a CVC relation between electromagnetic and weak form factors 
in the isospin limit. However, both the aforementioned theoretical accuracy and the new 
high-precision experiment at Brookhaven warrant a closer investigation of isospin 
violation. A crucial quantity in this connection is the pion mass difference M 2 + — M 2 = 
0.067 M 2 that is almost exclusively due to electromagnetic effects. Therefore, both the 
light quark mass difference and electromagnetism have to be taken into account in a 
consistent treatment of isospin violation. 

We concentrate in this letter on isospin violation in the reactions t~ — > tt°7t~v t 
and e + e~ — > 7r + 7r~ at low energies. Chiral perturbation theory (CHPT) || is the only 
framework where such corrections can be reliably calculated for the standard model in a 
systematic low-energy expansion. More specifically, we are going to calculate the leading 
corrections of both 0[{m u — m d )p 2 } and 0(e 2 p 2 ) for the CVC relation between the two- 
pion (vector) form factors in the two processes. A systematic chiral counting will be 
essential to extract the leading effects at low energies. 



2. The contribution of hadronic vacuum polarization at 0(a 2 ) to the anomalous mag- 
netic moment of the muon = (g^ — 2)/2 is given by || 

< CP ° l = ^ J 4M2 dtK(t)a° e+e ^ dmns (t) (1) 

where K(t) is a smooth kernel concentrated at low energies. The superscript in 
c"e+ e - ^hadrons denotes the "pure" hadronic cross section with QED corrections removed 
For the two-pion final state under discussion this means that Fy(i) in 

*° e+ e-^At) = na2P l + t At) \Fy(t)\ 2 (2) 

/W(t) = A 1/2 (l, M 2 + /t, Ml + /t) = (1 - AMl+lt) 1 ' 2 

X(x,y,z) = x 2 + y 2 + z 2 — 2(xy + yz + zx) (3) 

is the vector form factor of the pion with QED turned off (except for electromagnetic 
contributions to the charged meson masses). 

The decay r~ — > 7r°7r~z/ r is in general governed by two form factors / + ,/_. In the 
absence of electromagnetic corrections, these form factors are functions of the single 
variable t that is again the invariant mass squared of the two pions in the final state. The 
inclusion of electromagnetic effects generates shifts to the form factors which depend on 
a second Dalitz variable u = (P T — p w ~) 2 . Denoting by f+(t), f-(t) the u-independent 
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components of the form factors (to be defined precisely below), the decay distribution 
with respect to t takes the general form 
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The factor Sew takes into account the dominant short-distance electroweak corrections 
||. In the discussion of semi-leptonic r decays, the QED scale of Sew is usually chosen 
at the t mass. Thus, to lowest order in a, the short- distance enhancement factor is 
given by Sew = 1 + ( a /^)k3g(Mf /m^). Including the dominant electromagnetic higher- 
order effects, one finds the commonly used value || Sew = 1-0194. The factor GW(£) 
arises from the integration of the w-dependent electromagnetic correction over the Dalitz 
variable u. In principle, the spectrum distortion Gem(£) receives both virtual and real 
photon contributions. The measured electronic decay rate of the r lepton is related to 
rf ) by § 
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In the isospin limitQ we have M w + = M w o, Sew = Gem(^) = 1 and 
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implying the CVC relation 
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Including isospin violation to leading order, 0[(m u — m^p 2 ] and 0(e 2 p 2 ), we find from 
Eq. (f|) that still only the form factor f+(t) survives to this order. The modified CVC 
relation takes the form 
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Bremsstrahlung of soft photons (in principle contained in the function Gem(*)) is 
subtracted (at least in some approximation to be discussed below) directly from the raw 
data ||. In the analysis of Ref. |J, some additional isospin- violating corrections such as 
the width difference T p + — T p o were applied. The importance of the phase space correction 
factor /3%+ % - (t) / /3%o n - (t) has very recently been emphasized in Ref. |ID[ . 

It is the purpose of this work to estimate the remaining contributions to Rmit). 
Working at leading order, the form factor Fy(t) needs to be calculated to 0[{m u — md)p 2 ] 
with physical meson masses (but without explicit photonic corrections) whereas /+(£) 
must be calculated to both 0[{m u — m d )p 2 } and 0(e 2 p 2 ) if dT/dt is to be extracted from 
actual t decay data. 



3. To first order in isospin breaking and to first non-trivial order in the low-energy 
expansion, isospin violation manifests itself in the pion vector form factor Fy(t) only in 
the masses of the particles contained in the loop amplitude: 



F v (t) = 1 + 2H n+n -(t) + H K+K -{t) 



(10) 



with 11 



H PQ (t) 



1 

F2 



hp Q (t,ii) + -tLl(fj) 



where F denotes the pion decay constant in the chiral limit. The expression for the loop 
function h P Q(t,fi) is reported in the Appendix. The low-energy constant Lg(fJL) governs 
the charge radius of the pion which is in turn completely dominated by the p resonance. 
We will use the prescription of Ref. |T2| where the CHPT form factor ( PD of 0(p 4 ) was 
matched to the resonance region: 
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with the hadronic off-shell width (for the present case of the p°, the charged pion and 
kaon masses must be inserted) 
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and with a subtracted loop function (setting \i = M p ) 

Reh PQ (t,M p 
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F 2 



(13) 
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The representation ([T2] ) has the following attractive features ||12|| : 



By construction, it has the correct low-energy behaviour to 0(p 4 ) and its asymptotic 
behaviour is in accordance with QCD; 

It gives an excellent description of e + e~ — > 7r + 7r~ data up to t ~ 1 GeV 2 with the 
single parameter M p ~ 775 MeV. 
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Figure 1: Photon loop diagrams (without wave function renormalization). 

For our present purposes, the representation ( |12|) exhibits in addition the correct be- 
haviour to first order in isospin violation. To the order we are working, neither the p + — p° 
mass difference nor isospin-violating corrections to F n (we use F = F n = 92.4 MeV) en- 
ter. All the isospin violation to this order is contained in the physical meson masses in 
T p {t) and HpQ{t) and this feature will carry over to the form factor f+(t), except for 
additional purely electromagnetic corrections. 

At our level of precision, p — u mixing does not appear either. Such higher-order 
effects are not necessarily negligible numerically (see, e.g., Ref. Jl3]]). They can be and 
partly are taken into account in the actual analysis of the data (see, e.g., Ref. 0) or can 
be included in the theoretical error of ay, acpo1 . 

4. To first order in isospin violation, this time including explicit photonic corrections, 
the form factor f + (t,u) is given by 

/+(*, u) = l + 2H^- (t) + H K o K - (t) + f?X(u, M 7 ) + . (15) 

Compared to the form factor Fy(t) in QTOP, the appropriate meson masses appear in the 
loop amplitude and there is an additional electromagnetic amplitude containing both the 
photon loop diagrams shown in Fig. [l] and an associated local part. The electromagnetic 
amplitude depends on the second Dalitz variable u = (P T — p n -) 2 but not on t. Using a 
small photon mass M 7 as infrared regulator, the electromagnetic amplitudes are given by 
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(the corresponding calculation for K i3 will be presented in Ref. 
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with more details) 
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Expressions for the functions A(u), B(u), and C(u, M 7 ) are given in the Appendix. We 
have included logarithmic terms in /j^j to make it scale independent. The coupling 
constants iT 12 , Xi, and X 6 appear in the low-energy expansion of the standard model at 



order e p with inclusion of virtual photons [JT5J and of both virtual photons and leptons 



T6fl . Here we have pulled out the short- distance part Xf D of X 6 by the decomposition 



13 



x r M = xi» + x r M 
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The size of these contributions is discussed in the next section. 

The loop function f\^,(u, M 7 ) encodes universal physics related to the Coulomb in- 
teraction between the r and the charged pion. In other words, the w-dependence of the 
loop amplitude has little to do with low-energy QCD and thus with the chiral expansion. 
Rather, it represents the contribution to the loop integral given by low-energy virtual 
photons (the ultraviolet part being absorbed in the definition of the local amplitude). It 
is therefore natural to factorize these universal effects in an overall term fT7| . Moreover, 
since this factorization does not rely on chiral counting, we are lead to write (cf. Eq. (|T2|)): 
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As in the case of Fy(t) in (0), this representation of f+(t) has the correct low-energy 
behaviour to 0(p 4 ) and it interpolates smoothly to the resonance region. The resonance 
width T p (t) in (|T3"D has to be calculated now with the appropriate h~tt° and K~K° 
thresholds and phase space factors. 

The photon loop amplitude f\o^(u, M 7 ) is infrared divergent depending on an artificial 
photon mass M 7 . This dependence is canceled by bremsstrahlung of soft photons making 
the decay distribution in (t, u) infrared finite. The sum of real and virtual contributions 
produces the following correction factor to the (t, u) decay distribution 



A(t, M ) = l + 2/ 1 c 1 >,M 7 



+ fi f brems(^,M,M 7 ,E 7 n 



(22) 
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which depends on the minimal photon energy E™ m detected in the apparatus and is 
independent of M T This factor has to be multiplied by the kinematical density D(t,u) 
(defined in the Appendix) and integrated over the variable u to produce the term Gem(^) 
in the decay distribution (f|) with respect to t: 

/■Mmax(t) 

/ duD(t,u)A(t,u) 

G ™(t) = '-^f^ • (23) 

/ duD(t,u) 

"'«min(t) 

The details of soft photon emission (and the function g hrems (t,u, M^, E™ m )) depend 
on the specific experimental setup. To the best of our knowledge, all r decay experiments 
relevant here || apply bremsstrahlung corrections in the same (approximate) way de- 
scribed in Ref. ||18|| : only the leading term in the Low expansion (proportional to 1/E 1 
in the amplitude) is taken into account including also the logarithmic term in the loop 
amplitude M 7 ) depending on M 7 (contained in the function C(u, M 7 ) given in 

the Appendix). As emphasized in Ref. ||18|| , this is only an approximate treatment of 
bremsstrahlung that can be trusted for sufficiently small E™ m . Assuming this prescrip- 
tion, the setup-independent part of A(t,u) therefore involves only the subtracted loop 
amplitude (x is defined in the Appendix) 

f^M = &u, M 7 ) + £ { Ml + ml -u) log, log J^L ( 24 ) 



replacing C p (u,M 7 ) in 

Following Ref. |17| , we have factored out a universal loop amplitude Zio^subW- Al- 
though this factorization is independent of the low-energy expansion it is interesting to 
analyse the dependence on the lepton mass m T . Unlike in K[ 3 decays [p4| , the charged 



lepton is not light compared to a typical hadronic scale ~ M p . However, we can perform 
an expansion in p/m T in complete analogy to heavy baryon CHPT where p stands for a 
typical meson mass or momentum. Expanding /i™ su ij(«) in inverse powers of m T yields 

£U<«) = £(- 1+1< *5§ + ° ( ^) • (25) 

It turns out that the function Gem(^) is quite insensitive to whether it is calculated from 
the full fx^ps^u) in ([24] ) or from its large-m T approximation fl2"5|) . The difference is 
negligible in the full range 4M 2 < t < 0.8 GeV 2 : the leading term in (^) provides an 
excellent approximation. 

5. The results of our analysis are summarized in Figs. |^ {a),{b) where we plot the 
function -Rib(^) and its component factors defined in Eq. (^) for 0.2 < t < 0.8 GeV 2 . We 
note that the dominant contribution at low t is given by the kinematical term @\_/0q_ 
| TO]. Photonic corrections embodied in Gem(^) reduce Rmit) in addition by about half a 
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Figure 2: (a) Correction factor -Rib(^) for isospin violation. The bands around the central 
curve correspond to the uncertainty in the low-energy constants (solid lines) and in the 
bremsstrahlung factor (dashed lines). (6) The separate factors defining Rm(t) in Eq.(Q) 
are plotted as solid line for [3\_/(3l_ ) dashed line for \Fy(t)/ f + (t)\ 2 and dotted line for 



percent, largely independently of t. The form factor ratio \Fy(i) / f + (i)\ 2 is dominated by 
the width difference T p + — T p o. 

We have used the following input for the calculation of Rmit). 

• We employ the form factor Fy(t) given in Eq. (|12|) . 

• We use the form factor /+(£) as given in Eq. (|2T|). The local contribution depends 
on three low-energy constants appearing in the chiral expansion. For the constant 
K[ 2 (n) a sum rule representation is available [H5|. Saturating the sum rule with 



low-lying resonances and choosing the QED renormalization scale between 0.5 and 
1 GeV, we arrive at the following estimate: 

K[ 2 {M P ) = -(3±1) x 1(T 3 . 

As for X\ and Xq(/i), no estimates are presently available. We therefore use the 
upper bounds suggested by dimensional analysis: 

\Xi\<j^, \X r e (M p )\< 



(An) 2 ' 1 bV - (4tt) 2 ' 

In the case of Xg(M p ), we have enlarged the naive estimate by a factor of 5 (the (3 
function associated with Xq ||16|| ). This accounts for the present uncertainty in the 



S 



matching to the short-distance contribution to X 6 performed in Eqs. ([TJ, [L9|). The 
corresponding uncertainty in Rm(t) is shown in Fig. ^| (a) (solid curves). 



We include the factor Gem(^) according to the discussion following Eq. (|22|). Due 
to neglect of sub-leading terms fl8[ in the function g^ Tems (t,u, M^, E^ m ), GW(t) 



can receive extra contributions of order a/(4n) x 0(1). We therefore assign an 
uncertainty of ±a/n to it. The effect on Rmit) is also shown in Fig. |](a) (dashed 
curves). Clearly, better knowledge of the radiative amplitude can be used to improve 
our determination of Gem(^)- 

In order to quantify the impact of Rm{t) on a™ cpo1 , we construct the following ratio: 



'^dtK(t)a^Z^(t)R lB (t) 

nt^) = , (26) 

where cjg+^,Y c +7r+7r _ (t) is obtained via CVC from the r decay distribution as given in 
Eq. (|7|). A few representative values of 7£(t max ) are given in Table To translate the 
ratio 7£(t max ) into a modification of a™ cpol (4M 2 <t< t max ), we take for the purpose of 
illustration the value a™ cpol (4M 2 < t < 0.8 GeV 2 ) = (4794.6 ± 60.7) x lO' 11 of Ref. [|. 
Let us assume for simplicity that this value is calculated from r decay data only. As 
the quoted number contains the ratio of the radiatively corrected hadronic rate and the 
measured electronic mode (including the radiative channel with a photon in the final 
state), we have to multiply it ]2(J by the correction factor 1 + (a/27r)(25/4 — n 2 ) of 
Eq. ([D in addition to 1Z(t max = 0.8 GeV 2 ). In this way, isospin violation would reduce 
vac P oi( 4M 2 < t < o.8 GeV 2 ) by 76 x lO" 11 or about one standard deviation of the 



a 



reported error. 



Table 1: Correction factor for a™ cpo1 due to isospin violation (defined in Eq. (|2~6D ) for 
some values of t max - An uncertainty of 0.002 - reflecting the one in the bremsstrahlung 
factor Gem(^) - should be assigned to the values reported here. This is also an upper 
bound for the uncertainty due to the low-energy constants (see Fig. 0(a)). 



tmax (GeV 2 ) 


0.3 


0.5 


0.8 


7^- (tmax) 


0.949 


0.974 


0.988 



6. We have calculated the leading isosp in-breaking corrections to the CVC relation 
between the e + e _ — > 7r + 7r~ cross section and the decay distribution for t~ — > 7r°7r~z/ T . The 
calculation was performed in the framework of CHPT to 0[(m u — md)p 2 )} and 0(e 2 p 2 ). 



9 



Our main result comes in the form of a function -Rib(£) displayed in Fig. §(a) that corrects 
the CVC relation for isospin violation and electromagnetic effects. Since -Rib(^) is smaller 
than unity in most of the region under consideration (4M 2 < t < 0.8 GeV 2 ) isospin 
violation accounts at least for a sizable part of the systematic difference at low energies 
between e + e~ and r decay data (e.g., Ref. 

In general, isospin- violating corrections are expected to be of the order 



A* 
Ml 



2 x 1(T 



(27) 



where M p stands for a typical hadronic scale and A„. = M 2 + — M 2 . Electromagnetic 
corrections embodied in the function Gem(^) are precisely of this magnitude as shown in 
Fig. 0(6) . In the form factor ratio \Fy(t) / f + (t) | 2 , the ratio ( p7| ) is enhanced by a numerical 
factor: 



F v (M i 



r 2 



1 + 



6A^ 



M 2 - 4M 2 



The biggest effect occurs in the phase space ratio ffl 



1.015 



(28) 



. It is governed by the function 



3A„ 



t-AM 2 



(29) 



and therefore dominates at low energies. 

Although the calculation is based on a low-energy effective description of the standard 
model we claim that the main features of the correction factor -Rfb(£) are valid up to 
t ~ 0.8 GeV 2 . Of the three factors in the definition (^j of Rmit), both the dominant 
phase space correction factor |I| and the photon loop effects are independent of the 
low-energy expansion. Finally, the main part of isospin violation in the form factor ratio 
\F v {t) / f + {t)\ 2 occurs in the p-width difference Y p + —T p o and should therefore be reliable 
in the vicinity of the resonance. 

In this work we have been concerned with corrections to the CVC relation between 
r data and the bare cross section (T® + + . An important next-to-leading-order effect 



of 0(a 3 ) in a^ ad involves final state radiative corrections in a e 



_^ 7T+n -. Within scalar 

QED, this was already calculated in Ref. |22[ and recently reported in [[HJ. At 0(e 2 p 2 ) in 
CHPT one obtains the same result |23] as in scalar QED because the local counterterm 
contributions cancel due to gauge invariance. The resulting correction to aj^ ad of 0(a 3 ) is 
positive 
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Appendix The loop function h P Q(t,fi) is given by 
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In terms of the variables 
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and of the dilogarithm 
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the functions contributing to f^,(u, M 7 ) are given by 
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The kinematical weight to be used in Eq. ( 



is 
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T -{m 2 T -t)+ 2M 2 M 2 - - 2 u (m 2 - t + M 2 + M 2 _) + 2 u 2 . 



D(t,u) 

The integration limits in Eq. ( PHj ) are given by 

Mmax/min(t) = o(t) ± b(t) 



a(t) 
bit) 



mz 



-A, + 1 + 2M 2 _+m 2 -t 
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